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Abstract: Extra-dimensional models can involve bulk pseudo-Goldstone bosons (pGBs) 
whose shift symmetry is explicitly broken only by physics localized on branes. Reliable 
calculation of their low-energy potential is often difficult because it requires an understanding 
of the dynamics that stabilizes the geometry of the extra dimensions. Rugby ball solutions 
provide simple examples of extra-dimensional configurations for which two compact extra 
dimensions are stabilized in the presence of only positive-tension brane sources. The effects 
of brane back-reaction can be computed explicitly for these systems, allowing the calculation 
of the shape of the low-energy pGB potential, V4,d{'{>)-, as a function of the perturbing brane 
properties, as well as the response of both the extra dimensional and on-brane geometries to 
this stabilization. If the (/^-dependence is a small part of the total brane tension a very general 
analysis is possible, permitting an exploration of how the system responds to frustration when 
the two branes disagree on what the proper scalar vacuum should be. We show how the low- 
energy potential is given by the sum of brane tensions (in agreement with common lore) when 
only the brane tensions couple to (p. We also show how a direct brane coupling to the flux 
stabilizing the extra dimensions corrects this result in a way that does not simply amount 
to the contribution of the flux to the brane tensions. The mass of the low-energy pseudo- 
Goldstone mode is of order ^ {fi/ F)'^mjfjf (where is the energy scale associated with 
the brane symmetry breaking and F < Mp is the extra-dimensional axion decay constant). 
In principle this can be larger or smaller than the Kaluza-Klein scale, rriKK, but when it 
is larger axion properties cannot be computed purely within a 4D approximation (as they 
usually are). We briefly describe several potential applications, including a brane realization 
of 'natural inflation,' and a dynamical mechanism for suppressing the couplings of the pGB 
to matter localized on the branes. Since the scalar can be light enough to be relevant to 
precision tests of gravity (in a technically natural way) this mechanism can be relevant to 
evading phenomenological bounds. 
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1. Introduction and Summary 

Brane-world models — where known particles are localized on surfaces within extra dimen- 
sions — have proven to be fruitful places to seek novel kinds of low-energy physics. As studies 
over the past decade show, their low-energy physics can be novel (relative to 4D models, say) 
because of several different mechanisms: 

Brane vs Bulk Kinematics: Because not all particles are trapped on the branes, different 
species experience the kinematics of different dimensions. This observation is what allows the 
existence of unusually large extra dimensions and low gravity scales [|l|, . 

Brane Back-reaction: Because branes can be localized within the extra dimensions, which 
need not be homogeneous, physical properties can vary from place to place within the extra 
dimensions. Such position dependence is generic once the back-reaction of localized branes 
onto their geometry is included 
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Dimensional stabilization: Because low-energy degrees of freedom can be lighter than the 
Kaluza-Klein (KK) scale, predictions for low-energy dynamics generically require an equally 
complete understanding of whatever physics stabilizes the extra dimensions |Q . 

These last two points in particular considerably complicate discussions of cosmology within a 
brane- world framework What appear to be shallow directions for the scalar potential on 
a brane for fixed bulk geometry can turn into much steeper directions once the bulk geometry 
is allowed to move. 

With a few exceptions, the interplay between brane back-reaction and the physics sta- 
bilizing the bulk remains relatively poorly explored. The main exception is the case of 
codimension-1 branes moving within one extra dimension, as for Randall-Sundrum (RS) 
models Q. Yet one wonders how representative codimension-1 systems are of the more 
generic situation having higher codimension, for which much less is known about brane back- 
reaction. Although some results exist for the back-reaction of branes on ten-dimensional 
geometries in string theory |^], it is often not possible to be as explicit about the form of the 
extra-dimensional geometry and its detailed interplay with brane back-reaction. 

In this paper we compute the low-energy potential for a class of codimension-2 brane 
models for which explicit compactifications involving both bulk stabilization and brane back- 
reaction are known. We focus in particular on the simplest of brane/flux compactifica- 
tions: (nonsupersymmetric) 6D Einstein-Maxwell theory with the extra dimensions stabilized 
through a competition between a background Maxwell flux and a bulk cosmological constant 
[0, |9|, To this system we couple a bulk Goldstone boson (axion), <f>, whose shift 

symmetry, (j) — )• 0-|-(constant), is explicitly broken by its couplings to the two branes that 
source the bulk through interactions of the schematic form 

Sb = J^(T,iO + <^b*J') , (1.1) 

where the integration is over the 4D brane world-sheet, E^, whose volume form is denoted to so 
Th{(p) represents a i;^-dependent brane tension. J-mn denotes the 6D Maxwell field strength — 
for which is the 6D Hodge dual — and the i;^-dependent coupling ^bi4') can be interpreted 
as the amount of Maxwell flux carried by the brane in question. 

By computing the back-reaction of the branes onto the bulk fields we obtain the low- 
energy potential for the resulting would-be Goldstone zero mode, if, that becomes a pseudo- 
Goldstone boson (pGB) ||T^ in the low-energy 4D effective theory. The branes back-react onto 
the bulk fields by changing their boundary conditions, through a codimension-2 generalization 



[|[^ , 14, 15| of the more familiar codimension-1 Israel junction conditions |16|. 

The resulting bulk field equations subject to the brane boundary conditions can be solved 
explicitly in some generality if the axion-dependence of the brane tension is regarded as a small 
change to a background, axion-independent value. In this limit the would-be zero mode is 
stabilized to a fixed value, <p = which we compute in two separate ways: first by explicitly 
solving the linearized field equations of the full 6D theory; and second by minimizing the 
dimensionally reduced axion potential in the 4D low-energy effective theory. Both methods 
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agree, and the generality of our result allows us to follow how the stabilized value and its 
energy density vary as a function of the axion couplings to the two branes. In particular we 
see what happens when the branes differ in the value for the axion that they prefer. As a 
by-product we also compute how the geometry of the extra dimensions changes due to the 
presence of the axion-brane couplings. 

The calculation reveals the following generic features 

1. In the absence of brane fluxes — i.e. = — the low-energy 4D potential is very 
generally simply given by the sum of tensions, summed over the branes present, 

v,s(.^) = Y,n[Mv)], (1.2) 

6 

where (pbiv) denotes the value taken by the (suitably renormalized) 6D scalar field at 
the corresponding brane position, regarded as a function of the zero mode ip. This 
agrees with the probe-brane approximation (which ignores brane back-reaction) since 
for brane tensions the contribution of the back-reaction first arises at second order. In 
particular the stabilized value, (/? = c/j^, satisfies 

Because the quantity {dTh/ d(t))^p^ governs the coupling of the lightest mode, ip, to matter 
localized on the brane, these couplings tend naturally to turn themselves off for small 
fluctuations of (p about its ground state. This could provide a phenomenologically useful 
mechanism for naturally decoupling light bulk scalars from brane matter, along the lines 



of similar earlier proposals |17|. 



2. By contrast, nonzero brane fluxes contribute at linear order in two ways, that are similar 
in size. The flrst arises because nonzero background fluxes are required to stabilize the 
bulk geometry, *T = —Quj (with, as before, uj the volume form). Because of this the 
^b*^ term modifies the value of the brane action, giving an 'effective tension' 

Tb{^)=n{^)-QM^)- (1-4) 

The second contribution arises because quantization of total flux requires the amount 
of bulk flux to change in response to the presence of flux on the brane, leading to an 
additional energy cost over and above that measured by the difference Tb — Tb. Although 
the complete expression for Vcs{ip) that results involves an integration of with respect 
to (/?, the predictions for (p^ and g^s ■= Vcsi'P*) turn out to be relatively simple. The 
prediction for tp^ is again given by eq. (|1.3| ) — with no contribution from — while 
the prediction for q^q becomes 

Qee = ^{TblMv^.)] - Q'^blM^.)]} = ^{nlMv.)] - 2Q<fb[M¥^.)]} . (1.5) 
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We see in this way that the brane flux 'contributes twice' to the vacuum energy at low 
energies. In particular, g^^g = for branes satisfying Tb[(j)b{(pi,)] = 2Q and so 

3. Assuming a 6D kinetic energy of the form F'^d'^' cf) dM4>i in order of magnitude the mass 
of the light 4D would-be zero mode, predicted by the low-energy potential is of order 
~ [fi/ F)'^mKKf{^i,), where ttikk denotes the KK mass scale and the (/^-dependent 
part of the brane tension is assumed to be of order In most cases of interest 
is given by 5^Tb/9</>^ evaluated at ip^,. When fi <^ F this mode satisfies <C rriKK, 
allowing its properties to be described in the effective 4D theory. (Unlike for a purely 
4D theory it makes sense to call (p a pseudo-Goldstone field even if ^ > i*", since its 
shift symmetry is everywhere unbroken in the extra dimensions except at the positions 
of the branes. However, if /u > F then the mass of the would-be zero mode becomes 
comparable with other KK modes, precluding calculating its properties within the 4D 
theory.^) 

The low-energy 4D theory obtained from models with bulk axions generically includes a 
light scalar field, whose small mass is technically natural because of the weakly broken shift 
symmetry. What the brane construction potentially provides is a UV completion that can 
explain why the masses and couplings are small in the first place. This could prove useful for 
a variety of low-energy applications, such as to extra-dimensional inflationary models, some of 
which we briefly describe below while examining specific examples of our general expressions. 

We organize our detailed discussion as follows. Our main results are presented in the 
next section, which starts in §2.1 by setting out the field equations describing the system of 
interest. These are then solved by finding solutions that are perturbatively close to simple, 
well-known rugby-ball solutions involving two branes interacting with a spherical 2D bulk. 
§2.2 does so first for the simpler case where the branes do not couple to the bulk scalar, 
with the generalization to scalar-brane coupling following in §2.3. §3 then explores the fea- 
tures of these general solutions by examining in detail several simple illustrative special cases. 
These include situations where the two branes agree on the value at which the field ip sta- 
bilizes, and situations where they do not. §3 also considers several special cases of potential 
phenomenological interest for axion and inflationary applications. 

2. The bulk-brane system 

This section defines the system of interest, which we take to be the simplest theory containing 
both codimension-two sources (with positive tension) and a flux mechanism for stabilizing the 
size of the extra dimensions. This suggests taking the bulk theory to be 6D Einstein-Maxwell 
gravity coupled to the Goldstone boson (axion) field. 

^Exceptions to this can arise if the scalar is self-locahzed at the brane but this usuaUy requires a bulk 
scalar potential U{(j>) that is forbidden in the current examples by the assumed shift symmetry, tf) ^ (p + c. 
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2.1 Field equations and background solutions 

Classical brane-bulk dynamics is defined by solving the bulk field equations, subject to the 
boundary conditions imposed by matching conditions at each brane. 

Bulk field equations 

The bulk action of interest is^ 

1 / 1 _ 1 . , „ A 1 



5buik = -y d^xV^ ^^g'"" i^-^n,,j, + ^J,,cl)d^<pj + - TmnJ''''' + , (2.1) 

where TZ^n is the Ricci tensor constructed from the 6D metric gMN, <i> is the axion field and 
T = dA is the field strength for the Maxwell potential Am- The dimensionful parameters 
of the problem are the 6D gravitational coupling, k = 1/Mg, the bulk axion decay constant, 
Ka = and the bulk cosmological constant, A (whose value is tuned to ensure the 

unperturbed solution is flat in the on-brane directions). Although Ka can be absorbed into 
the normalization of (j), we do not do so because this changes the form of the brane couplings 
to (p. 

The fleld equations obtained from this action are the (trace-reversed) Einstein equation 



TZmn + a dM4>dN(t> + k FmpJ^n^ 
where A := k/kq = F/Mg. The Maxwell equation is 



Kj A 



gMiv = 0, (2.2) 



^F~9 Vm-F*^^ = dM [y^g ^^^^ j = , (2.3) 

and 

^gU4> = Bm W^g d"cj)) = , (2.4) 

is the axion equation. 
Rugby-ball solutions 

We consider geometries that are maximally symmetric in the 4 on-brane directions and axially 
symmetric in the two extra dimensions. The corresponding ansatz for the metric, scalar and 
Maxwell fields is 

ds^ = dp2 + e'^^dO^ + e^^ gf.udx^dx" , (2.5) 

and 

J'pe = A'e , (2.6) 

where g^u is an x'^-dependent maximally symmetric geometry and the functions B, W, (/) and 
Ag depend only on p. Primes denote differentiation with respect to this coordinate. 



We use a 'mostly plus' metric and Weinberg's curvature conventions (that differ from those of MTW 
pc| ] only by an overall sign in the definition of the Riemann tensor). 
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Subject to this ansatz the bulk field equations reduce to 
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(H, (2.7) 



where 72. is the curvature scalar built from the maximally symmetric metric g/j^i^. The first 
two of these immediately integrate to give 

0' = (^1 and 6"^+^^^^ = Q , (2.8) 

where ipi and Q are integration constants. 

When Lfi = the full set of equations admits a particularly simple solution of the rugby- 
ball form 

ds^ = dp2 + a^L^ sin^ d^^ + da^'^dx'" 

Tpe = QaL sin , (2.9) 

with constant (p = (po and W = 0. The equations of motion imply the following relation 
amongst the integration constants: 

2 ^2(^^A^ (2.10) 



L2 V 2 
as well as fixing the 4D curvature 

n := t'T^^u = (Q' - 2A) . (2.11) 

A final constraint relating parameters comes from flux quantization, due to the spherical 
topology of the extra dimensions. As usually framed, this implies 

^ = 2aL'^Q, (2.12) 

where g is the gauge coupling of the Maxwell field and n is an arbitrary integer. However 
this expression assumes the absence of any flux localized on the source branes themselves 
[|lC| ] . In the presence of brane-localized flux (more about this below and in Appendix ^) the 
flux-quantization condition instead becomes 

^ = 2aL^Q + Y.P', (2.13) 
g ^-^ 27r 



- 6 - 



where the sum is over all of the branes present, each of which carries the localized flux, 

Since our interest is in background solutions with flat geometries, g^y = rj^y, we further 
choose A so that = 0: 

A = ^ and so k^L^Q^ = 1 . (2.14) 

With this choice all geometrical properties, like L and Q, can be regarded as functions of the 
integration constant a together with the integer n and the lagrangian parameters k and g 
(see Appendix ^ for details). 

The potential singularity in the geometry where gee = vanishes is just a coordinate 
artefact when a = 1, in which case the two compact dimensions define a sphere. When 
a 7^ 1 the background has a conical singularity at p = p^f := and p = ps '■= ttL. This 
is interpreted as describing the back-reaction of two codimension-two source branes located 
at these positions, having equal tensions, T (which includes the energy associated with the 
brane flux, see eqn. ( 2.16| ). Matching at the branes (see below) implies this tension is related 
to the deficit angle by 

l-a = ^, (2.15) 



similar to the relation between tension and deficit angle for a cosmic string [21|. 

Finally, notice that the value </) = 930 is not determined by any of the equations of motion, 
due to the symmetry (/>—)•(/> + constant. The parameter ^pQ labels a flat direction in the low- 
energy potential, that can be lifted if the coupling of (j) to the branes breaks this symmetry 
(such as by allowing the tensions T to depend on 0). A primary goal of the next few sections 
is to identify the effective potential for this low-lying mode below the KK scale, to determine 
how the vacuum value after symmetry breaking, c/p*, is related to the couplings on the branes. 



Brane matching conditions 



As brane sources we use the most general form (involving the fewest derivatives) for a 4D 
brane action located at positions p^ and ps |llO| 



'branes 



b=N,S 



' b 2 '^'^ 



b=N,S 



mn 77 

'6 n ^ mn 



(2.16) 



where e^^ = Xj^fg^ = transforms as a tensor, rather than a tensor density, in the two 
transverse dimensions. The parameter t^, represents the tension of the brane, which can 
depend on all of (j), W and gee without breaking the condition of maximal symmetry in the 
on-brane directions. As is shown below, the parameter similarly denotes the magnetic 
charge (or flux) carried by the source branes (which could also depend on cp, W and gee)- 
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The presence of such branes imposes a set of boundary conditions on the derivatives of 
the bulk fields in the near-brane limits, given by'^ 



d% 



with Th 



and 



1 



Pb 



Pb 



Ub with Ub :- 



Tb + 3Ub 



-2 ^dTb 



4vr \dgee 



(2.17) 



where Tb is defined as the total lagrangian density of the source, 

Tb = Tb- ^be-^Fpe. 



(2.18) 



The Bianchi identities ensure that only two of eqs. ( p. 17 ) are independent of one another, 
and as a consequence the quantities Ub and % are also not independent. They are subject to 
the constraint: 



2 -2Tb- 3Ub 



(Tb 



/^2 



0, 



(2.19) 



where 7^' = dTb/dcj). Notice that for the rugby ball solutions =Ub = Q and so eqs. ( p. 17] ) 
degenerate down to eq. ( 2.15| ). 

As shown in Appendix the corresponding boundary condition for the Maxwell field 
implies that the integral of eq. ( |2.8D for Ae{p) for a patch containing each source brane is |1C] 



Mp) 



2tt 



Q 



dp Northern hemisphere 



Pn 



— + Q fdpe^-'^^ Southern hemisphere , (2.20) 
2vr 



and the signs are dictated by the observation that increasing p points away from (towards) 
the North (South) pole, together with the requirement that the two patches share the same 
orientation. Requiring these to differ by a gauge transformation, g~^dg^l, on regions of overlap 
implies the flux-quantization condition 



n 

g 



tot 



27r 



PS 



dp e 



B-iW 



(2.21) 



PAT 



where n is an integer, g is the gauge coupling and <I>tot = ^jv + ^5- It is this expression that 
identifies ^b as the fraction of the total Maxwell flux carried by each brane. 

2.2 Perturbations I: the Einstein-Maxvifell case 

Next consider starting with a rugby-ball solution and independently perturbing each of the 
two brane tensions, T{, = r + 6Tb, and brane-localized fluxes, = <^ -|- d^b^ implying a similar 



^Notice that we normalize the quantities % and Ut without inchiding a factor of e*^ used in ref. pA 
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expansion for the total brane action, Tf, = Tf, — e^J-po^h = Tb — Q^b- This section starts 
simply and assumes both 6ti, and are independent of (j), with the resulting insights used to 
inform the next section's discussion of the more general case. The goal is to compute explicitly 
how the bulk fields respond to the perturbation, allowing a detailed examination of how the 
extra dimensions flex as their source branes change. In general, because the perturbed branes 
are different from one another, the scalar field acquires a nontrivial profile, (j) = (j){p), and the 
resulting geometry warps nontrivially, W = W{p). 

Linearized solutions 

Because the brane perturbations are independent of cp, dTh/dcp = and so there is no change 
to the (j) boundary conditions. Consequently the unperturbed solution, = ipo, remains a 
solution. The scalar then drops out of the problem and the calculation involves only the 
Einstein-Maxwell system. Writing = e^°[l + (5i?(/?)] and W = 6W{p) — with Bq given by 
the rugby-ball solution described by parameters Q, L and a — we linearize the field equations 
in 5B and 6W. 

The combination of the Einstein equations (pp) — {00) linearizes to 

5W" - B'q 6W' = 6W" - ^ cot (^) = (2.22) 



L \L 
which has as its solution 

dW = Wi cos , (2.23) 

where we absorb an additive integration constant, Wq, into a re-scaling of the four on-brane 
coordinates, x'^. 

The perturbed gauge field again satisfies eq. ( p.8D , and so 

With this and eq. ( p.23| ) the 66 Einstein equation linearizes to 



5Q + Q{SB -45W) aLsin(£) . (2.24) 



whose integral is 



Notice that the integration constant Bi here is pure gauge, corresponding to an infinitesimal 
shift in the radial coordinate /? — )• p + c. We can fix this freedom by demanding that /Ojv = 0, 
and so — >• as p — )■ 0. Since = aL{l + 6B) sin{p/L) — t- —Bi at /? — )■ 0, this implies 
Bi = 0. 



The linearized flux quantization condition, eq. ( 2.21 ) is 



- 9 - 



which uses the background relation nLQ = 1 to rewrite 5^>f,/L^Q = K^Q5<^f,- Solving this for 
5Q/ Q gives 

-± = -iSBo ^ i6<Pr, + S<^>s) . (2.28) 

y na 

In summary, once coordinate conditions are used to eliminate Wq and Bi, solutions to 
the bulk equations for 6W, 6B and 6Ag involve three integration constants — Wi, 6Bq and 
SQ/Q — among which flux quantization imposes one relation. The physical interpretation 
of the two remaining parameters is seen by connecting them to two physical quantities. One 
of these can be taken as the warping difference between the two branes, 



6W^ - 6W. = 2Wx 



(2.29) 



which controls the relative redshift of energies on the two branes. The other can be chosen as 
the change in proper distance, ps — = vr(L + 5L), between the two branes, where /O^y = 
and ps are defined as the places where — 0. Comparing 



lim 



SnaL 



6Q 
Q 



with the Taylor expansion e^(7rL) ~ {e^")' {—tt6L) = +7raSL gives 



6L 
T 



6Q 
Q 



(2.30) 



(2.31) 



Matching to brane tensions 

All that remains is to eliminate the integration constants 6Bq and Wi in terms of the brane 
perturbations using the linearized brane matching conditions. In the present instance only the 
last of eqs. (|2.17] ) is nontrivial. Besides imposing the background relation a = 1 — K^T/(27r), 
for the linearized perturbations this condition implies 

Kt 27r ' 



6 e 



(2.32) 



where 6Tb = 5{n - ^bQe'^^") ^ 6n - 6<^bQ - '^Q{SQ/Q - 4.5Wb). Evaluating this at 
P = Pn = and^ p = ps = ttL + 5ps, and keeping in mind that it is —p that is the outward 
direction for the south brane, gives 



a 



a 



■ „ WWi 3 



6Q 
Q 



a 



a 



-25Bq - 
-25Bo + 



lOWi 


3k^Q6^ tot] 


K^5T^ 
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Ana 


271 


lOI^i 


3K2Q5^>tot' 


K^6Ts 
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Ana 


2tt 



(2.33) 



When solving these we may approximate 6Tb — 5Tb — Q6^b^ which involves dropping the 
back-reaction of those terms proportional to 5Q/Q and 5Wb in r?bTbl2'K. This neglect is 



*To leading order ips does not contribute, and is only mentioned for completeness. 
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justified because their relative contribution is of order K^Q<I>/27r, which must be small to 
justify our classical treatment of gravity. The solution found within this approximation to 
eqs. ( 2.33 ) then is 



Svra 
3^2 



6Tj, + 5Ts - 3Q((5$^ + (5$c 



<5r„ - 5T, 



(2.34) 



Using the above, 



6L 
T 



6Q 
Q 



3k^ 
Svra 
3^ 
Svra 



(5r^ + 6Ts 



{5tn + 5ts 



2Q(5$^ + 6<^s] 



(2.35) 



On-brane geometry and the view from 4D 



The curvature of the induced geometry on the branes comes from the linearized {fj.i') Einstein 



equation, which 



n 



using e '^^ F^g 



-8W 



5W" + 



cot 



L 



- gives 



+ 



L2 ' L2 



Q 



2 

l2 



5_Q 
Q 



(2.36) 



vraL^ 

From the point of view of a 4D observer localized on the brane this curvature would be 
interpreted as being due to a 4D energy density, £>eff- Since the 4D gravitational coupling. 



SvrGjv, is related to k by 



1 27raL 



we have 



Ka 



n 

'44 



ttL 



dp sin ( — 



bTj, + 5ts - 2Q(5$jv + 



(2.37) 



(2.38) 



Notice that this agrees with the naive expectation ^)eff = '^''"jv + brg in the absence of fluxes 
on the brane. The same is not true in the presence of brane fluxes, however, since the final 
result for ^>efT differs from bT^^ + 5Ts = br^ + Sts — Q{b^N + b^s)- As the above calculation 
shows, ^eff —bQ/Q and so the energy cost of the perturbation arises from the change of 
flux required by the flux-quantization condition in response to the back-reaction of the branes 
on the bulk geometry. Since the flux is homogeneous across the extra dimensions, its energy 
cost is expensive since it scales with the volume. Localizing some of the flux into the branes 
reduces this extensive energy cost.^ 



^Of course, this possibility of back-reaction competing with brane tensions is aheady suggested by the 
complete absence of on-brane curvature in the initial rugby ball solution despite the presence of the initial 
equal brane tensions, T. 
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The comparative importance of such back-reaction effects depends on the relative size 
of the two brane energy scales Sti, and Q6^b = d^t/nL. If both 6ti, and 6^h are set by the 
same scale — i.e. Tb ~ and S^b ~ A — then 6Tb ~ QS^b when A ~ A^. := (kL)^^^^. 
For A smaller than this the Q5^b term dominates, while 5Tb is the larger of the two when 
A > A*. (For a similar setup with n transverse dimensions this crossover would occur when 
A ~ (kL)"^/^^"''"').) Although kA^ must be much smaller than one to justify semiclassical 
methods, for fundamental objects it is comparatively large {e.g. for D-branes kA^ is of order 
the string coupling, gg ~ 0.01 say), and so the tension contribution can therefore dominate. 
The flux contribution instead can dominate for lower-tension objects. 

It is instructive to check this calculation by directly evaluating the low-energy potential 
through dimensional reduction of the 6D theory in the classical approximation. A general 



formula for this is computed (including brane back-reaction) in ref. |14], and when this is 
specialized to linear perturbations about a rugby ball it evaluates to 

V., = 2./"dpe— [?^c„t(|)] -l(Q + .Q)V-4} , (2.39) 

with the W term arising from the extra-dimensional curvature and the (Q + (5Q)^ term com- 
ing from the bulk Maxwell action. In the present instance all terms involving 6W in this 
expression turn out to be proportional to sin(p/L) cos(p/L) to linear order in the perturba- 
tions, and so integrate to zero and do not contribute to VeS- Keeping in mind the background 
relations 2 A = and kLQ = 1, the result therefore simplifies to 

showing the equivalence between the 4D and 6D perspectives. 
2.3 Perturbations II: the Einstein-Maxwell-axion case 

In this section we generalize the previous discussion to consider branes and fluxes that can 
depend on (p. This allows us to follow how couplings to the brane lift the flat direction 
associated with the shift symmetry of the bulk theory, and so to see how the scalar zero 
mode, ifQ, becomes stabilized at a specific value, (/Jq = V'*- 

It is instructive to ask how this stabilization happens from the point of view of the full 
six-dimensional theory. To this end imagine trying to integrate the field equations to obtain 
the bulk configuration that interpolates between the two branes. Specializing to solutions that 
are both axially symmetric in the transverse directions and maximally symmetric in the on- 
brane dimensions we seek bulk profiles as a function only of p, starting with initial conditions 
set by matching to the brane at p = (say). If this matching completely specified all of the 
fields and their first derivatives at this brane then the solution obtained by integration would 
completely determine the value of the fields and their radial derivatives at the second brane, 
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and in general these need not be consistent with what would be obtained by matching to this 
second brane. 

But matching to the first brane typically specifies only the derivatives of the fields at the 
first brane, and not separately the values of the fields themselves.^ Consequently the values 
of the fields at the first brane can be adjusted to try to allow the solution to properly match 
to the properties of the second brane. It is in this way that the system can force 999 = if 
the brane actions do not preserve the bulk shift symmetry. 

From the perspective of a low-energy 4D observer the energy cost responsible for this 
stabilization looks like a scalar potential for (/jq, and our goal in what follows is to compute its 
shape for configurations in the immediate neighborhood of 999 = 92*- As the above arguments 
show, a classical solution subject to our assumed ansatz should not exist as soon as (/pq 7^ V*) 
and the part of the ansatz responsible is likely to be the condition of maximal symmetry in 
the on-brane directions. No maximally symmetric solution should exist for LfQ 7^ 93^, because 
this indicates the onset of time evolution in response to no longer sitting at the minimum of 
the 4D effective potential. (This development of time dependence in response to changes in 
the properties of mutually gravitating brane sources resembles what happens for a system of 
electric charges, which generically becomes time dependent when an equilibrium arrangement 
is disturbed). 

Rather than trying to solve for the system's time-dependent response (see however 
refs. [22, ^]) when 939 7^ 99* we instead focus on computing features of the low-energy poten- 



tial that is responsible. We do so - in both the 4D and 6D theories - through the artifice of 
turning on a current that stops the time evolution, and so removes the obstruction to static 
solutions. Since sufficiently small deviations from equilibrium should precipitate motion only 
along the low-energy fiat directions, it suffices to couple the current only to these low-energy 
modes. By computing the amount of current required as a function of the low-energy scalar 
mode, we may Legendre transform in the usual way to determine the shape of the effective 
potential. 

Linearized equations with currents 

Since we work within a linearized approximation, we perturb the brane properties in a way 
that does not drive the low-energy scalar fields far from their initial values. This can be 
achieved if the potential energy of each brane has a minimum as a function of </>, and although 
the two branes need not agree on where this minimum is they should not disagree by too much. 
It suffices therefore to study the brane tensions in the vicinity of these minima, restricting to 
quadratic expansions in powers of 4>. Writing Tf, = T + 5Th{<p), we take 

5n{^) = no + ^{(i)-vt)\ (2.41) 

with b = N and S and Ti,(</)) = r{,(</>) — Q^bi^P), as before. For technical reasons — see 
Appendix |D| — we require that the minimum of the sum of the brane actions, = 0, 

^Since in general the bulk fields can diverge at the brane positions, this argument should more precisely be 
made very near to, and not precisely at, the position of the first brane. 
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agrees with the minimum of the sum of the fluxes, Ylb ^'b ~ When ^ Vs the two branes 
differ on which value for (p they prefer, and we assume that this difference is not so large as 
to invalidate a linearized integration of the field equations. 

In the higher-dimensional theory the current used to stabilize the solutions against rolling 



is^ 



d^x V=5 J{<P - V>^) , (2.42) 



and, to the extent that it suffices to stabilize just the KK zero mode, J can be taken to be 
independent of the extra-dimensional coordinates p and 9. In the presence of such a current 
the field equation for (j) becomes 



Om [V^ g'^'d^ct)) =^KiJ. (2.43) 

Perturbing around the rugby ball solution, our interest is in the lowest nontrivial order in 
J, corresponding to situations where the brane tensions only cause controllably small changes 
in (p. In this case the leading approximation to the axion fluctuation is obtained by solving 
eq. ( |2.43| ) with the metric evaluated at the rugby ball background. 



sin(^)50' =llsin(^), (2.44) 



L 



e.7 . fP 



L2 \L 



where the last equality defines the dimensionless current, assumed small: ej := k^JL^ <^ 1. 
The Maxwell equation is unchanged by the current, and integrates to give 

Fpe = A', = {Q + 5Q)e-4^+^ . (2.45) 

Recall that both W and B in this expression include perturbations. 

Including the stress energy from the current interaction, the linearized Einstein equations 
become 

,sw" , SB" . ^ cot (I) . ^H^'f , -L (f ) - 5^ , . „ 

- ^ - (!) - - (f ) - ^ (f ) - ^ - ^ » . (-) 

where, as before, X = k/ Ka = kF. 



^The additional coupling J ip* is here inserted to ensure that J couples only to the light scalar mode, 
Sip — ipo ~ ifii,, at the linearized level, and not also to the metric fluctuations. We keep this term even though, 
as discussed in Appendix for axions much lighter than the KK scale, m <^ a misalignment that 

included metric modes only introduces subdominant contributions to the axion mass, of order Sm^ ~ rn^L^ . 
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Linearized solutions 



We now solve those equations to order €j. Because {'p')'^ is order e^, to order ej the equation 
for the warping is unchanged from previous sections, giving the solution 



6W = Wi cos 



(2.47) 



The current forces the axion to acquire a profile (which is desirable because this allows it 
to satisfy the new boundary conditions at the brane positions). The perturbed axion equation 
integrates to give 



S(t) = (po + ipi In 



1 — cos(p/L) 



sm{p/L) 



e,,ln 



(2.48) 



with ipi and (po integration constants. Since ipo parameterizes the (previously) flat direction 
we solve €j and all other integration constants in terms of it and brane properties. 
Using these in the Einstein equations as before gives SB as the solution to 



[SB'sin^ip/L)]' _ lOWi 



sin^ (p/L) 



6Q 
2L2 V Q 



2L2 



(2.49) 



givmg 



6B 



3 fSQ\^Xhj, 



^cot(|)-l^cos(|)...o. (2.50) 



Using this in the linearized flux-quantization condition finally gives a relation between 
SQ and SBq, 



SO k'^O 

= Xhjiifo - if^) - 4SBo (5$^ + 5$s) 



(2.51) 



As before, the remaining integration constants — in this case (pi, W\ and SBq — are deter- 
mined by solving the matching conditions at the brane positions. The fractional change in 
the proper distance between the source branes becomes 



SL_ 

T 



3 (SQ\ X^ej 



4 V q) 



(2.52) 



Matching to branes 

The matching condition for the axion at each brane is 



lim ap 4) = 
p-^Pb 27r 



(2.53) 



p->pb 



but an additional complication arises because the right-hand side is ill defined due to the 
divergence in (j){p) at the brane positions. This requires a renormalization of the parameters 
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defining the brane potentials []l5| , 24]. To this end first regularize the matching condition by 
evaluating it at p = + and p = pg — £3- Then define the renormalized parameters 



Vn ■■= Vn - ifi - ej) In (^^^ + ln2 

991 In 2 . 



Vs ■= Vs + ivi + ej) In ( j- 



e . In 



Because the field profile satisfies 

(and a similar result at p = vrL — eg), these definitions ensure 

4){eb) -Vb = ^o-Vb, 



L 



(2.54) 



(2.55) 



(2.56) 



and so remain finite in the limit — t- 0. This makes the derivative of the tension (and the 
tension itself) finite when evaluated on the brane. Because the fluxes are also written in terms 
of — Vh, they do not need a separate renormalization. 

In terms of renormalized quantities the matching conditions directly relate the integration 
constants, 



-a{<pi + €j 



27r 



Kl5T'g{ipo) (kITs2 



2tt 



2tt 



2tt 



{^0 - Vn) 
{^0 - Vs) , 



(2.57) 



allowing the inference 
,^2 



Ana - 



5T'f,{ipQ) + 5T's{ipQ) 



Aira 



{Tn2 + Ts2) v?o - Tn2 Vn - 7s2 Vs 



Ana - 



Aira - 



{Tn2 - Ts2) (fO + Ts2 Vs - Tn2 Vf, 



(2.58) 



The flrst of these identifles the held value where the fiat direction gets stabilized, ipQ = (p^,, 
since this is the solution that corresponds to zero external current. The condition ej{ipi,) = 
implies satisfies 

5T'^[^,)+5T'g{ip,) = Q, (2.59) 



and so when 5Tb{ipQ) = + i Tb2 (y'o - Vb)"^ 



Tn2Vn + Ts2Vs 



Tn2 + Ts2 

We again fix 6Bq and Wi from the last of the matching conditions, eqs. ( |2.17] ), 

1 



V J n, 



T + 6Tb{ipo) 



(2.60) 



(2.61) 
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which uses Ub{(pQ) ~ 0, as can be inferred either from the second of eqs. ( p. 17 ), or by solving 
eq. (2.1£) to hnear order in k^T^. As before this leads to the conditions 



a 



a 



-25 Bn 



-25Bo + 



tot 



3 47ra 
WWi 3K;2Q5$t 



The result is 



Aira 



3k^ 



j.2\ 



^^j6Ts{vo)- (2.62) 



(2.63) 



and 



\ 2 ^ 

SBo = ^(V^o - V.) + ^{ST^i^o) + STsiipo) - 3Q [5^^(990) + S<^s{^o)] } (2.64) 



.2 
Svra 



^ I X] ['^^b('^o) - (V'o - ^*)6Tl,{^o) 

\^b=N,S 



3Q5^tot(v'o) 



where the final line eliminates Ej using eq. ( ^.581) and A^k^ = k^. This implies 

J.2 

27ra 

^2 



2 

-XhA^o - V.) - ^{'^^^(V'o) + STsiifo) - Q5^tot{vo)} (2.65) 



27ra I ^ 

b=N,S 



1 



STb{ipo) - -(v?o - 93*) '57fe(¥'o) 



In terms of these the fractional change in the proper distance between branes becomes 
- = -41 



L 



2 
8na 



3k^ r ^ 



Q6'^toti<fo) 



(2.66) 



b=N,S 

On-brane geometry and 4D effective potential 

The linearized Einstein equation yields the following on-brane curvature 

86W 2 /6Q\ 2X^tj 



7^ 



2 

I2 



m" + ^ cot [P- 

— - A ej{(po - (p^. 



12 ^ L"^ \Q 



L2 



(2.67) 



ira 



L2 



b=N,S 
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The presence of the current J compUcates the determmation of the effective potential, 
T4ff ((^), in the low-energy 4D theory. The appropriate matching calculation turns on a current 
in the low-energy theory as well, and asks what potential reproduces the previous results for 
TZ and 

The most general action for the 4D effective theory involving only the 4D metric, g^u^ 
and the low-energy scalar, is (up to two derivatives) 



5. 



eff 



d X V-^ 



1 
24 



(2.68) 



where j is the low-energy current, K4 is given by eq. ( p. 37 ) and the 4D axion decay constant 
is 

. . A-wnT? X2 

(2.69) 



= A-KaL^F* = — = — . 



We couple the current j to the difference 93 — (/j^ purely as a matter of later convenience. 

The equations of motion, specialized to constant scalar fields, (p = (^Oi and to maximally 
symmetric geometries, are 



J 



-Kff(93o) and ^ = -j((/?o - P*) - Kfr(v5o) , 



from which j can be eliminated to give 

((^0 - ¥'*)Kff(V5o) - 14fr(v'o) 



4^1 



(2.70) 



(2.71) 



The functional form for the potential Ves is determined by requiring eq. ( 2.71| ) to repro- 
duce the curvature, eq. ( ^.67 ), predicted by the 6D theory, regarded as a function of tpo. This 
can be obtained by regarding eq. ( p. 71 ) as a differential equation for Vcs, whose solution is 



dip 



{if - ip^f 



(2.72) 



1^0 



dip 



Q5<^toti<p) 
if - f*f 



lim 



In general, the coefficient of the term linear in (po — ip^,) in T4g is the integration constant, 
which is fixed in the second equality of eq. ( p. 72 ) by requiring V^g((/?j,) = 0, for ip^ as given 
in the 6D theory by eq. ( p. 591 ). 

Two physical parameters of particular interest here are: {i) the effective on-brane cos- 
mological constant, ^>efr := V^ff(¥'*), and the low-energy scalar mass, := V^g((/9^)//^. The 
first of these evaluates to 



= 6Ti^{ip^) + 6Ts{ip*) - 2Q 6^N{pi,) + 6^s{y^* 



(2.73) 
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whose value agrees with our ear her Einstein-Maxweh calculation when JT^ = 0. Similarly' 



m,, 



P 



1 

7^ 



5T';,{ip,) + 5T'^{^^)+ lim 



(2.74) 



In the absence of brane fluxes the effective potential is simply the sum of the brane potentials. 
But although the low-energy scalar always stabilizes at the stationary points of 'Ylib^'^b^ the 
scalar masses and 4D cosmological constant in general differ from what would be expected 
based just on 5Tf,. 

Comparison with dimensional reduction 

As before, we can also evaluate V^s at the classical level by direct dimensional reduction, 
which gives the integral 









'8W' 








\2k2 


[ L 





ttL 







dp . fP 

T Kl 



2TTa (5Q ,2 



TvaL 



Q 



n. 



(2.75) 



in agreement with the 6D calculation above. 



3. Applications and special cases 

This section seeks to illustrate the physical implications of the previous section's results by 
exploring several instructive examples. 

3.1 Bulk response to stabilizing potentials 

Consider first the response of the bulk geometry and the properties of the low-energy 4D 
scalar-tensor theory, distinguishing the cases where the two brane agree on, or compete for, 
the field value where the low-energy scalar is stabilized. 

Shared minima 

As an example where the fluxes and tensions on both branes are minimized at a common value 
of ifo, consider the special case that all the fluxes and tensions have the following expansion 
T^ = T + 5n{<j)) and = $ + 5^ij{(j)) with 

5n{ipQ) = 6TM + ^^{ipo-vf and 5^h{po) = ^^bo + ^^^{po - vf , (3.1) 

^Notice that diverges if $tot vanishes hnearly with ip — ipi,. In this case the lowest energy KK mode is 
not properly captured by our ansatz — see Appendix |^ — and so the low-energy potential misidentifies its 
size. It is for this reason that we require $tot and Ttot to agree on the value ipi, at which they are minimized. 
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and so 



5T, 



(3.2) 



where ST^k = Srhk — Q S^bk are constants. 

The condition fixing ip^, in this case is ip^ = v, as one would expect. Inserting this into the 
formulae for the relative warping of the two branes and the fractional change in inter-brane 
distance gives 

3/^2 



6Wn - SWs 
SL 

T 



207ra 



STno — dTso 



(3.3) 



Similarly, the on-brane expressions for ^eff aiid yield 



^eff = STt,o + SSTso - q(s^ MO + S^so) 



(3.4) 



and 



= ^ [^1^2 + STs2 + Q(S^n2 + 6^32) 

6tn2 + ^rs2 



f 

1 r 
P 



(3.5) 



Notice that only the second derivative of the tension, 5T'^{(pi,), contributes to the scalar mass, 
while both the tension, 6Tb{ipi,), and the flux, 6^b{ipi,), contribute to the on- and off-brane 
curvatures. 

Brane competition 

Consider next the case where the two branes each prefer (pQ to stabilize at different values, 
causing them to compete in the value they ultimately determine. A representative example 
in this case is 



ST. 



(5$ 



S^b{ipo) = 5^bo + {(fo - Vb 



(3.6) 



The stabilizing value for the scalar is now neither nor Vs, but instead the intermediate 
value 

5Tm2Vn + STs2Vs 
01 n2 + 01 s2 

with the ratio 5X^2 /STs2 controlling precisely where (p^, lies between and Vg- Requiring 
<^^tot = + S^s also to have its minimum at the same value of then requires 



S2 



ST. 



S2 



(3.8) 
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The extra-dimensional geometry satisfies 



5W^ - 5W. 



207ra 



5T, 



NO 



6Tff2STg2 



5Ts2 - ST] 



n2 



{6T^2 + 6Ts2y 



and SL/L = (Sk^ Qcs/Sna), with the 4D vacuum energy given by 



QeS 



1 / ST^2STs2 \ 



2 [6T^2 + STs2j 



Q 
2 



{6T^2 + STs2)^ 



{Vn - VsY 



1 



6Tn2 + STs2 



(3.9) 



(Vr.-Vs)'', (3.10) 



where the last equality uses eq. (p^). The result for is again given by eq. (|3.5|) . Because 
9?^ does not minimize the tension at either brane both the total tension and total flux get 
increased by positive amounts. These positive contributions then act oppositely in ^eflf- 

More complicated competitions can also occur if there is also symmetry-breaking in the 
bulk, in which case competition between the bulk and brane potentials can lead to self- 
localization [fql. 



Flux domination 

A particular instance of the previous scenario corresponds to the case where \5ti,\ <C |Q(J<J>fe|, 
since in this case STh{(p) — —Q6^b{(p). Then the stabilizing value for the scalar becomes 



and 



5W, 



' 20TTa 



5^ 



NO 



5^so + 



S^n2Vn + S^.S2V.S 
6^n2 + S^S2 ' 

6<^n2S^S2 \ ^^S2 



6^ 



n2 



{S<^n2+S<^s2V 



(f jv - VsY 



(3.11) 



(3.12) 



while SL/L = {SK^gcs/^T^Oi), with 
QeS '- 



{vm - VsY 



(3.13) 



Because in this case ~ to leading order, the scalar mass — and so also the stability of 
the vacuum (/?o = V'* — is controlled by subdominant effects (like 6Th or loops), even though 
the flux dominates the classical contribution to £»eff- 

3.2 Axions 

It is instructive to consider the relative sizes of the various scales that arise naturally when 
bulk axions receive masses through their couplings to branes, since these need not be related 
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in the same way as when both axion and symmetry-breaking physics share the same number 
of dimensions. This section briefly examines several illustrative choices. 

There are five scales that naturally arise in bulk-axion models. Three of these — the 
extra-dimensional Planck scale, Mg = k~^I'^\ the axion decay constant, F = Ka^^"^; and the 
KK scale, rriKK = l/L — characterize the bulk physics. The source branes are responsible for 
the other two: the scale A set by the (/>-independent parts of the brane tensions and fluxes; 
and the scale set by the (?!)-dependent terms, 

~ A\ $60 =i a , rfe2 =i / and $62 - A* ■ (3.14) 

These scales are not completely arbitrary. In general, control over the semiclassical 
approximation requires Mg to be much bigger than all of the others. Although the conditions 
Hih? <^ 1, Kfi^ <^ 1 and k/L"^ <^ 1 follow fairly directly from standard arguments [^], the 
condition F <^ Mg is a bit more indirect. Because F sets the scale of the bulk symmetry 
breaking for which <j) is the would-be Goldstone boson, our upper bound on F assumes the 
UV completion describing this breaking intercedes below the Planck scale (before which the 
UV completion associated with gravity — such as string theory — should also intercede). 

Furthermore, we generically expect /i A for generic types of brane physics. This follows 
because it is difficult to have physics contribute to the (p mass without also contributing 
equivalently to the vacuum energy. Notice in this regard that it is technically natural to 
take /i <C A, because it is only fi that breaks the shift symmetry of the low-energy scalar: 
if ^ if + (constant). 

Axion mass 

In terms of these scales the mass of the light scalar in the effective 4D theory is of order 

m^:^^=.(|)'i:^(|)'m^,, (3.15) 

in all three of the scenarios considered above. ^ This result doesn't depend on which scenario 
is considered because for all three the scalar mass depends only on \/6Tpf2 + ^Ts2/f- Provided 
fi <^ F, if is much lighter than the KK scale as is appropriate for its description in the 
low-energy 4D eff'ective theory. 

For the higher-dimensional models of interest here, however, the regime fj,^ F can also 
make sense. Extra dimensions allow this regime even though the scale /i of explicit symmetry 
breaking is then much larger than the scale of the spontaneous breaking: F. Because all 
symmetry breaking is localized on the branes, even though ^ > F the field (/> behaves like 
a Goldstone boson for all energies lower than F in the bulk provided one stays away from 
the position of the branes. Although this regime is not amenable to a 4D description, the 
mass of all KK modes can be computed within the higher-dimensional theory. In this limit 

®In some circumstances additional suppression can be achieved, such as if the Goldstone symmetry is not 
completely broken by either brane separately 
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the 'zero mode' becomes lost among the generic massive KK states and is not singled out 
as being particularly light. In this regime it is clear that otherwise standard arguments, 
like cosmological bounds on axion properties, cannot be made purely within four dimensions 
without taking the full dynamics of the extra dimensions into account. 

Curvatures 

A second robust prediction of all of the above scenarios is the relation between the change to 
the extra-dimensional size and the four-dimensional curvature: 

although the size of g^s itself is not as model independent. This source of this model de- 
pendence is the competition between tension and flux contributions to g^s, whose competing 
contributions are of order j£»eff — 'l2b'^bo or dgcs — ^5 Q^bo, with 

A M^A , (Mf,A\ 

no ~ A^ and Q^feo ~ QA ~ — ~ ~ — ^ . (3.17) 

kL L \ Mp I 

Special things happen for the BPS-like situation when the tension and charge are precisely 
related, r;,((/?*) = 2Q^f,{ipi,), since in this case the two contributions to g^s precisely cancel. 

Whether the tension or the flux dominates in g^g depends on where A sits relative to 
the two geometrical scales Mp ~ 10"'^'^ GeV and 1/L. Defining A^ := yf^jwiMg /Mp we have 
^^eff — A^ if A > A^ and ^^cff — AA^ when A < A^. Some representative numerical values are 
given in Table 1. Intriguingly, A^ is of order the QCD scale in the extreme case of large extra 
dimensions {Mg < 10 TeV and m^A- < 0.4 eV ||). 



Mg 


A. 




10^5 


4 X 10^^ 


4 X 10^2 


10" 


2 X 10^ 


4 X 10^ 


10^ 


8 X 10^ 


4 X 10"^ 


10^ 


0.8 


4 X 10-1° 



Table 1: Values oi tukk and A^ as a function of Mg (in GeV). 
3.3 Gravitationally coupled scalars 

The special case F ~ Mg is of particular interest because then / ~ Mp and the low-energy 
4D scalar is gravitationally coupled. In this case the light scalar mass is robustly of order 
— jj?' /Mp, and its small size is technically natural since it is protected by the underlying 
shift symmetry. There are two situations for which the existence of such light weakly-coupled 
scalars are of particular interest. 
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An inflationary mechanism 

Inflationary models famously require light, weakly coupled scalars; something that is usually 
fairly difficult to achieve without fine-tuning in a real microscopic theory. The above estimates 
point to a fairly generic mechanism for achieving slow-roll inflation whenever a bulk axion 
acquires a potential through its interaction with codimension-two branes. This mechanism 
can be regarded as an ultraviolet completion of 4D 'natural inflation' models that assume 
the inflaton to be a pseudo-Goldstone particle. 

The mechanism rests on two assumptions: (i) the brane energy density, Qeg, must domi- 
nate any other contributions to the geometry in the on-brane directions; and (ii) the brane- 
axion couplings must have a local maximum rather than a minimum at (p = if^,, for which 
is of order fi^ /Mp (as above) but negative. In this case because the previous estimates apply 
near the potential's maximum, with an effective 4D scalar potential being of order 

Veeiip) A + BU (if - if^) , (3.18) 

with B ~ 0{fj,'^) and A = g^s ^ 0{A^) or g^s ~ 0{AM^/Mp) (whichever is larger). The 
technically natural choice /u ^ A, (Mg /Mp)^^^ therefore ensures B <^ A. Here U{x) is a 
calculable, dimensionless, order-unity function, whose expansion for small arguments is (by 
assumption) U{x) ~ — ^ ?72 + • • • with U2 > and order unity. 

Should this potential dominate the 4D geometry it produces a Hubble scale near this 
maximum that is of order H ~ y^^g/Mp and so H is of order the larger of /Mp or 
{M'^/Mp){K/MpY'^. Because of this, our choice fi <^ A automatically ensures |m^| <^ H. 
Provided that H is also small compared with the KK scale — as is easy to arrange — 
the resulting cosmology can be understood within the 4D effective theory, and describes an 
inflationary slow roll provided ip starts in an initially spatially homogeneous configuration 
near the potential's maximum. This slow roll is inflationary (despite having / ~ Mp) because 
B <^ A, since the slow-roll parameters are of order e ~ {BU' /A)^ and 77 ~ BU" /A. rj is 
sufficiently small to inflate for ~ 60 e-foldings ii B/A 0.01, in which case e ~ r/'^ is even 
smaller (and so the inflation typically does not produce an observable signal of primordial 
gravity- waves) . If ~ then B/A ~ (/u/A)^ and a sufficiently small ratio can be ensured 
for the comparably modest hierarchy fi/A ~ 0.3. 

As an existence proof that all parameters can be chosen as required above consider the 
intriguing, but extreme, scenario where the QCD axion is a bulk scalar within large extra 
dimensions (the last line of Table 1). In this case taking fi ~ A/3 ~ Aqcd — 0.2 GeV both 
provides the right scale of axion-matter couplings, and ensures A ~ A^ and so ges ~ A"^ 
and H <^ ninK- One might imagine that whatever solves the cosmological constant problem 
arranges the true ground state of the present epoch to be the unperturbed rugby-ball solution 
having 7^ ~ and T of order the weak scale, with the perturbation Jr^ ~ Ag^^ arising in the 
early universe due to the vacuum energy associated with the QCD phase transition on the 
brane. Even if it were not to involve enough e-foldings to account for primordial fluctuations, 
such a very late inflationary period could be useful for removing unwanted relics — like moduli 
or KK modes — from the much earlier universe. 
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New long-range forces 

Another potential application (or constraint) on the light bulk Goldstone mode described 
here comes from the long-range forces that it would mediate if its mass is sufficiently light. 
Indeed, one motivation to study the brane-bulk dynamics explored above is to find sensible 
UV completions which can have a technically light scalar whose presence could be sought 
when testing general relativity. Such tests provide strong constraints on the existence of any 
new forces competing with gravity in the solar system, with a precision that varies with the 



mass of the new scalar particle and the nature of its couplings to matter [28|. 

This section explores what brane-bulk dynamics might say about the couplings of the 
low-energy scalar to matter localized on the branes. We find these couplings can (but need 
not, depending on the brane properties) realize some earlier-proposed mechanisms for 
dynamically vanishing when the scalar is in its ground state. 

To see how ip couples to matter localized on the branes, we generalize the previous discus- 
sion to include brane- localized matter fields, generically denoted by ifj. All brane quantities 
like tension and fiux are regarded as being functions of both brane and bulk fields, 

n = nitp, (t>), ^b = ^b{ip, (t>) and so on. (3.19) 

The main point is that none of this afi^ects the matching conditions and solutions described 
above, and so in a static (or adiabatic) configuration the ground-state value = v^*(^) still 
adjusts to satisfy 



E 



= 0. (3.20) 



The new ingredient that appears in searches for new forces is the use of spatially inho- 
mogeneous matter configurations as sources {e.g. planets, stars, etc.) of spatial variation, 
Sip = 5ip{x), for the fiuctuation 5(p = ip — (p^, along the on-brane directions. Regarded graph- 
ically, these constrain the amplitude for emitting a single (p particle from the source, with 
repeated emissions accumulating to give a coherent classical field. But the amplitude for 
(/^-emission from matter localized on a specific brane, h = b^, is controlled by the expansion 
of the brane action in powers of the fiuctuation. 



Tboii^,'^) = T6o(V',</?*) + 



(5(/? + C'(V). (3.21) 



Of these interactions, it is only the term linear in 6ip that acts as an obstruction to solving 
the field equations with 5(p = 0, and so it is this linear term that is subject to the strongest 
constraint from new-force searches. Unless T^q has special properties such a term could 
generate violations of the equivalence principle, which are strongly excluded once the range 
of the force becomes macroscopically large. 

Comparing eqs. ( |3.20| ) and ( |3.21 ) reveals the mechanism for suppressing (/9-matter cou- 



plings. If the action, T^g, for the specific brane on which we live should share the same 
extremum as does the sum of all branes, Th, then as adjusts to satisfy the condition 
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( |3.2[1| ), it would also automatically turn off the dangerous coupling of 5^p to matter localized 
at brane h^. As the examples above show, the extremum of the sum of all brane actions need 
not agree with the extrema of each brane's action separately. But it automatically does so in 
two simple cases: (z) when none of the branes besides ho couple to (f) at all; and {ii) when all 
of the branes couple to </>, but are all extremal for the same place. 

Notice that the argument is not changed by the presence of (/)-dependent brane fluxes, 
<I>fe. This is because they do not enter into eq. ( p. 20 ) independently from their contribution 
to Tfe (even though they do contribute independently to the value of ^'eff)- 



Acknowledgment s 

We wish to thank Markus Luty and Raman Sundrum for useful discussions about back- 
reaction in codimension-two models. CB acknowledges the Kavli Institute for Theoretical 
Physics in Santa Barbara and the Abdus Salam International Center for Theoretical Physics 
for providing the very pleasant environs in which some of this work was performed, as well 
as Eyjafjallajokull for helping to provide some unexpected but undivided research time. LvN 
thanks the Instituut-Lorentz for Theoretical Physics at Leiden University for their hospital- 
ity. Our research is supported in part by funds from the Natural Sciences and Engineering 
Research Council (NSERC) of Canada. Research at the Perimeter Institute is supported in 
part by the Government of Canada through Industry Canada, and by the Province of Ontario 
through the Ministry of Research and Information (MRI). 



A. Brane fluxes and flux quantization 

To see how to interpret the parameter rewrite the brane flux term as a regularized 6D 
integral weighted by a scalar function s{p) whose support is nonzero only in a short interval 
1/9 — Pb\ < £ away from the brane, and is normalized so that J d'^x s = 1. That is, 

5flux = Y jd^x 7^ s e^^:Fmn = j d^x s Fpe . (A.l) 
Then the SAe Maxwell equation becomes 

dp FP' - $b s) = , (A.2) 

which integrates to give 

e^^{e-''Mg-^ks) = Q. (A.3) 

This is the bulk solution found in the text away from the brane, where s = 0. 

Imagine now integrating this to obtain Aq{p) in the vicinity of the brane at = 0, using 
for s a simple step function: s = l/(7re^) for p < e and s = for p> e. Assuming W ~ Wi, 
is approximately constant and ^ p foi p < e, the solution satisfying .40(0) = is 
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and so at p = e in particular 



(A.5) 



The junction condition for A'g at p = e can also be seen by subtracting the solution, 
eq. ( [A.3| ) evaluated at p < e — where s = 1/ (vre^ ) — from the same solution evaluated at 
p > e — where s = 0. Since the RHS is the same in both cases we get the following jump 
discontinuity across p = e: 

e-^A'g =--^- (A.6) 
L J p=e- vre^ 

This can be related to the derivative of the brane action with respect to Ae by rewriting 
eq. ([Ol) as 

Sflux = ^6 y d^'x 7^ s F^e = ^ j d^x Ae{e) , (A.7) 
and so (keeping in mind the relative sign between the tension and flux terms) 



P=e+ 1 
p=e- 2tt V dAt) 



(A.8) 



as stated in ref. [pH] . 



B. Rugby-ball response 

This section provides the explicit solutions for the properties of rugby ball solutions as func- 
tions of the assumed (shared) brane tension, and in particular computes the response to small 
changes in its value. 

Rugby ball configurations solve the field equations 



-77 -pPQ 



9mn — . 



and 



subject to the ansatz 



ds^ = g^u dxf'dx" + dp'' + a^L^ sin^ 



L 



TpQ = aQLsin 



L 



(B.l) 
(B.2) 

(B.3) 
(B.4) 



with (j) = ipo constant. The bulk field equations give the 2D and 4D curvature scalars as 



<2) = ^ 



3Q2 



+ A 



and 



2lQ' 



7^(4) = n = 2K^^- A 



(B.5) 
(B.6) 
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The gauge potential corresponding to eq. (B^) is 
^± = aQL^ 



±1 



p 

cos ( - 



27r 



(B.7) 



where the it sign indicates the solution for the northern or southern hemisphere, since A must 
evaluate to the brane flux at the corresponding pole. Requiring the difference between these 
two solutions near the equator to be a well-defined gauge transformation, gA^ — gA~ = dO, 
implies the constants Q and L must be related by 



gQ 



N 



(Bi 



where we define N = n — (7$tot/2vr. 

Eqs. ([B.5[) and ([B.8[) determine the constants Q and L in terms of a and A, with solutions 



1 



2^2 



Sa^g 



L2 3iV2^2 



and 



1 ± 



3iV2^4A 

8 a^g"^ 



1 



1 ± 




2a5L2 ~ 3ArK2 




(B. 



(B.IO) 



These provide two solutions for L and Q for each given value of a and A, satisfying L? > 
-^min ~ K^/ 16 a^g"^. Starting with the lower sign in eq. ( [B.9D the radius L falls from 
L — )• oo to L = \/2 Lmin as A climbs from to A^ax = Sa'^g'^/SN'^K'^. On this branch 
A< A 

max implies ~ k?A/2. Then switching to the branch corresponding to the upper 
sign has L fall from Lmin to Lmin as A recedes from A^ax back to zero. There are no real 
solutions with A > Amaxj or with L < Lmin- 

For each of these solutions the last equation, eq. (|B.6| ), gives the on-brane curvature, TZ. 
There is a choice A = Aj, for which it vanishes, given by Aj = Q^/2. For this choice L and 
Q become 

1 2ag 



Lf{a) Nk 



and Qf{a) 



2ag 



and so 



2„2 



(B.ll) 



(B.12) 



2 iV2^4 • 

Because Lmin < < -^min we see that this solution lies on the branch corresponding to 
the upper sign of eq. (|B.9|) . In particular 

2 / AT \ 2 

(B.13) 



2ag\' 



Nk 



2ag 



1. 



Notice that the semiclassical approximation requires the curvature to remain small com- 
pared with the relevant energy scales, and so in 6D requires 7i? to be much smaller than A 
or Q2^ Because TZ ~ k2A and this requires k^A and must both be much smaller 
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than unity. So for A ~ Q ~ a g /N k the semiclassical hmit imphes a g /N k <^ 1. This 
in turn ensures k/Lj ^ 1, showing that this value of Ues within the classical limit. 

The above expressions can be used to check the linearized analysis performed in the main 
text. To this end, suppose we start with a = uq, with A = Aq = Aj(ao) chosen so that TZ = Q 
for this value of a. Then we change the brane tension (but not the brane flux), and so also 
a, without also adjusting A. Choosing the upper sign, the radius and magnetic flux become 



1 

l2 



2^2 



Sa^'g 



1 + Wl 



3Af2K4Ac 



2a' 



1 



l + \ 1 



Sag 
4a2 



(B.14) 



and 



Q 



1 + 



802^2 



2a \ 
3ao/ 



1 Qo 



1 + A 1 



Sag 
4a2 



(B.15) 



The last equality in these two equations is obtained by using eq. (B.12) to trade Aq for ao, and 



then using eqs. ( B.ll ) to express the result in terms of the values Lq and Qo that correspond 
to a = ao. These equations show how the values of L and Q adjust to compensate for the 
change of a. The on-brane curvature similarly changes, and is given by 



9A^2^2 



2Ao) 



4a2 



2a2 



(B.16) 



which vanishes as a — )• ao, as it must. For a = ao + Aa, then ag/a ~ 1 — 2Aa/ao and so 



7^ 



(\^alg^\ /Aa\ _ 4 /Aa\ _ 
\W^) \ ~ Z2 \ ~ ~ TTa^Ll 



SkIAT , 



(B.17) 



which uses the matching condition 1 — a = 4GT = k^T/Itt in the form Aa = —k'^AT/2it, 
as well as the definition of the 4D gravitational coupling: k'^ = AttuoLq k^. Defining the 4D 
potential, Vcs, hy TZ = —^k\ Fcff gives the expected result 



Feff ^ 2Ar . 



(B.18) 



The factor of 2 arises because a change of a requires an equal change of tension for both 
branes if it is to preserve the rugby-ball form. 



C. Misaligned currents 

This Appendix uses a simple model to track the implications that arise if the external current 
happens not to be aligned precisely with the lightest mode of the system. When this happens 
errors can arise in the identification of quantities like low-energy masses, but this section 
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argues that these are generically suppressed by powers of the light mass divided by heavier 
masses. 

Consider then the toy 4D lagrangian 



-9 



'2^2^ 2 



whose potential is given by 



(C.l) 



(C.2) 



with masses assumed to satisfy m ^ M. Here 99 is meant as the analog of the KK would-be 
zero mode in the main text, while x is representative of some other, more massive, KK mode. 
The goal is to ascertain the extent to which our method of determining the low-energy mass 
would be thrown off by a small coupling — parameterized here by C — of the external current 
to a heavy state. 

The classical equations of motion for the scalar fields are 



□ - m^^ = J and □ x - M'^x = C J : 



while the Einstein equation reads 



'R'fiiy + Ki ( df,ip duip + d^x duX ) + = , 



and so 



7^ 

4^ 



Evaluated at the particular solutions 



J = -m If and x = -jp 



M2 



this last equation gives 



1 + 



(C.3) 

(C.4) 
(C.5) 

(C.6) 
(C.7) 



In terms of J-{^) the method of the main text gives the low-energy scalar potential as 

VM^):=ipl ^H^)- (C.8) 



For T{ip) = A + Bip + ^dp'^ the integral evaluates to 



10 



es = -A + Biplnip + - Cip^ + Dip , 



(C.9) 



^^The singular form of V(,'4(0) when B ^ corresponds to the pathological case where brane fluxes and 
tensions are not extremized for the same value of (fi, discussed in more detail in Appendix]^. 
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where D is the integration constant, and so when apphed to the above toy model this gives 

2" 



1+1^ 
' M 



(C.IO) 



This expression correctly identifies the value of the potential at its minimum to be Vq, and 
— provided ( < 0(1) — gives the correct mass for the field ip, up to corrections of relative 
order rn^/M'^. 

D. When brane fluxes and tensions compete 

This Appendix briefly discusses another kind of competition, which would arise if (^Tf,((/?o) and 
(5<I>fe((/9o) at the same brane were not minimized by the same scalar configuration. A simple 
representative in this category is 

T^{ipo)=T + T^o + ^iVo-VTf and <l>^i^o) = <^ + '^no + ^ i^o - v^)\ (D.l) 

together with 6Ts{(po) = 6^sifo) = 0, so the 'south' brane plays no role in the stabilization 
of (/Jo- 
in this case because the flux is irrelevant for determining ip^, its value is simply ip^ = Vt- 
The warping difference is also insensitive to $b and so becomes Wjs, — Ws = (3/t^Tjvo/207rQ), 
while 6L/L = {3k'^ Qcs/^iTa) with 

QeS = Tjvo - Q^jvo - -^y^ {vt - v^f 

and ™j = ^ + %i,im (^). (D.2) 

Clearly, the expression for the mass is singular when Vt v^. The reason for the singularity, 
is that for this choice of brane there is no solution satisfying the ansatz with which we work. 
The obstruction lies with the Maxwell field, which we choose to lie in the J-'pQ direction only. 
However, the perturbation that gets excited by moving ipo away from equilibrium, if we do 
not stabilize with an external current, necessarily gets a time dependent Maxwell field. But a 
changing magnetic field induces an electric field, so the Tpt and Tet components cannot both 
remain zero. 

To see that the Maxwell field must acquire a time dependence, consider a perturbation, 
5ip, that oscillates about the background Vt, 

ip = vr + difip) e"^""* . (D.3) 

In the fiux condition, eq. ( |2.51| ), the brane fiuxes now have a part that is linear in 6if that 
becomes proportional to e"*™*. If we now assume that 6Q does not acquire any time de- 
pendence, we get a contradiction: 6B has a part proportional to e"*™* according to the fiux 
condition, but in matching it with the brane in eq. ( |2.32| ), if 5Q is time independent the 
right hand side is either constant or proportional to — ^ e~^*"^*. This shows that the 
matching conditions cannot be satisfied unless the magnetic field becomes time-dependent. 



- 31 - 



References 

[1] N. Arkani-Hamccl, S. Dimopoulos and G. Dvali, Phys. Lett. B429 (1998) 263 
[hep-ph/9803315]; Phys. Rev. D59 (1999) 086004 [hep-ph/9807344]; 

I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B436 (1998) 257 
[hep-ph/9804398]. 

[2] K. Benakli, Phys. Rev. D60, 104002 (1999) [hep-ph/9809582]; C. P. Burgess, L. E. Ibanez and 

F. Qucvcdo, Phys. Lett. B447, 257 (1999) [hep-ph/9810535]; 

P. Horava and E. Wittcn, Nucl. Phys. B475 (1996) 94 [hep-th/9603142]; Nucl. Phys. B460 
(1996) 506 [hep-th/9510209]; 

E. Witten, Nucl. Phys. B471 (1996) 135 [hep-th/9602070]; 
J. Lykken, Phys. Rev. D54 (1996) 3693 [hep-th/9603133]. 

[3] L. Randall, R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370 [hep-pli/9905221]; Phys. Rev. Lett. 
83 (1999) 4690 [hep-th/9906064]. 

[4] W. D. Goldbcrgcr and M. B. Wise, "Modulus stabiUzation with bulk fields", Phys. Rev. Lett. 
83 (1999) 4922-4925 [arXiv:hep-ph/9907447] 

[5] P. Binetruy, C. Deffayet, D. Langlois, "Non-conventional cosmology from a brane-universe" 
Nucl. Phys. B 565 (2000) 269-287 [arXiv:hcp-th/9905012] 

P. Kanti, I. I. Kogan, K. A. Olive, M. Pospelov, "Cosmological 3-brane solutions" Phys. Lett. 
B 468 (1999) 31-39 [arXiv:hep-ph/9909481] 

J. M. Cline, C. Grojean, G. Servant, "Cosmological expansion in the presence of extra 
dimensions" Phys. Rev. Lett. 83 (1999) 4245 [arXiv:hep-ph/9906523] 

[6] B. R. Greene, A. D. Shapere, C. Vafa and S. T. Yau, "Stringy Cosmic Strings And 
Noncompact Calabi-Yau Manifolds," Nucl. Phys. B 337 (1990) 1; 

S. Kachru, R. Kallosh, A. D. Linde, J. M. Maldacena, L. P. McAllister and S. P. Trivedi, 
"Towards inflation in string theory," JCAP 0310, 013 (2003) [arXiv:hep-th/0308055]; 

D. Baumann, A. Dymarsky, L R. Klebanov, J. M. Maldacena, L. P. McAllister and 
A. Murugan, "On D3-brane potentials in compactifications with fluxes and wrapped 
D-branes," JHEP 0611 (2006) 031 [arXiv:hep-th/0607050]; 

D. Baumann, A. Dymarsky, I. R. Klebanov and L. McAllister, "Towards an Explicit Model of 
D-brane Inflation," JCAP 0801 (2008) 024 [arXiv:0706.0360 [hep-th]]. 

[7] A. Kehagias, "A conical tear drop as a vacuum-energy drain for the solution of the 
cosmological constant problem," Phys. Lett. B 600 (2004) 133 [arXiv:hep-th/0406025]; 

T. Kobayashi and M. Minamitsuji, "Brane cosmological solutions in six-dimensional warped 
flux compactifications," JCAP 0707 (2007) 016 [arXiv:0705.3500 [hep-th]]; 

D. Kiley "Rotating Black Holes on Codiniension-2 Branes," Phys. Rev. D 76 (2007) 126002 
[arXiv:0708.1016 [hep-th]]. 

[8] F. Leblond, R. C. Myers and D. J. Winters, "Consistency conditions for brane worlds in 
arbitrary dimensions," JHEP 0107 (2001) 031 [arXiv:hep-th/0106140]; 

S. M. Carroll and M. M. Guica, "Sidestepping the cosmological constant with football-shaped 
extra dimensions," [hep-th/0302067]; 



- 32 - 



I. Navarro, "Codimension two compactifications and the cosmological constant problem," 
JCAP 0309 (2003) 004 [hep-th/0302129]; 

E. Papantonopoiilos and A. Papazogloii, "Branc-bulk matter relation for a purely conical 
codimension-2 brane world," JCAP 0507 (2005) 004 [arXiv:liep-th/0501112]. 

[9] A. Salam and E. Sezgin, "Chiral Compactification On Minkowski X S**2 Of N=2 
Einstein-Maxwell Supcrgravity In Six-Dimensions," Phys. Lett. B 147 (1984) 47; 

[10] Y. Aghababaie, CP. Burgess, S. Parameswaran and F. Quevedo, Nucl. Phys. B680 (2004) 
389-414, [hep-th/0304256]. 

[11] S. L. Parameswaran, G. Tasinato and I. Zavala, "The 6D SuperSwirl," Nucl. Phys. B 737 
(2006) 49 [arXiv:hep-th/0509061]; 

H. M. Lee and C. Ludeling, "The general warped solution with conical brancs in 
six-dimensional supcrgravity," JHEP 0601 (2006) 062 [arXiv:hep-th/0510026]; 

H. M. Lee and A. Papazoglou, "Supersymmetric codimension-two branes in six-dimensional 
gauged supergravity," JHEP 0801 (2008) 008 [arXiv:0710.4319 [hep-th]]; 

C. P. Burgess, S. L. Parameswaran and L Zavala, "The Fate of Unstable Gauge Flux 
Compactifications," arXiv:0812.3902 [hep-th]; 

[12] C. P. Burgess, "Goldstonc and pscudo-Goldstonc bosons in nuclear, particle and 
condensed-matter physics," Phys. Rept. 330 (2000) 193 [arXiv:hep-th/9808176]. 

[13] C. P. Burgess, D. Hoover and G. Tasinato, "UV Caps and Modulus Stabilization for 6D 
Gauged Chiral Supergravity," JHEP 0709 (2007) 124 [arXiv:0705.3212 [hep-th]]; 

C. P. Burgess, D. Hoover, C. de Rham and G. Tasinato, "Effective Field Theories and 
Matching for Codimension-2 Branes," JHEP 0903 (2009) 124 [arXiv:0812.3820 [hep-th]]; 

C. P. Burgess, D. Hoover and G. Tasinato, "Technical Naturalness on a Codimension-2 
Brane," JHEP 0906 (2009) 014 [arXiv:0903.0402 [hep-th]]; 

[14] A. Bayntun, C. P. Burgess and L. van Nierop, "Codimension-2 Brane-Bulk Matching: 

Examples from Six and Ten Dimensions," New J. Phys. 12 (2010) 075015 [arXiv:0912.3039 
[hep-th]]. 

[15] P. Bostock, R. Gregory, L Navarro and J. Santiago, "Einstein gravity on the codimension 2 
brane?," Phys. Rev. Lett. 92 (2004) 221601 [arXiv:hep-th/0311074]; 

I. Navarro and J. Santiago, "Gravity on codimension 2 brane worlds," JHEP 0502 (2005) 007 
[arXiv:hep-th/0411250]; 

J. Vinet and J. M. Cline, "Codimension-two branes in six-dimensional supergravity and the 
cosmological constant problem," Phys. Rev. D 71 (2005) 064011 [hep-th/0501098]; 

M. Peloso, L. Sorbo and G. Tasinato, "Standard 4d gravity on a brane in six dimensional flux 
compactifications," Phys. Rev. D 73 (2006) 104025 [arXiv:hep-th/0603026]; 

B. Himmetoglu and M. Peloso, "Isolated Minkowski vacua, and stability analysis for an 
extended brane in the rugby ball," Nucl. Phys. B 773 (2007) 84 [hep-th/0612140]; 
B. Himmetoglu and M. Peloso, "Isolated Minkowski vacua, and stability analysis for an 
extended brane in the rugby ball," Nucl. Phys. B 773 (2007) 84 [hep-th/0612140]; 

E. Papantonopoulos. A. Papazoglou and V. Zamarias, "Regularization of conical singularities 
in warped six-dimensional compactifications," JHEP 0703 (2007) 002 [arXiv:hep-th/0611311]; 



- 33 - 



"Induced cosmology on a regularized brane in six-dimensional flux compactification," Nucl. 
Phys. B 797 (2008) 520 [arXiv:0707.1396 [hep-th]]; 

D. Yamauchi and M. Sasaki, "Brane World in Arbitrary Dimensions Without Zi Symmetry," 
Prog. Theor. Phys. 118 (2007) 245 [arXiv:0705.2443 [gr-qc]]; 

N. Kaloper and D. Kiley, "Charting the Landscape of Modified Gravity," JHEP 0705 (2007) 
045 [hep-th/0703190]; 

M. Minamitsuji and D. Langlois, "Cosmological evolution of regularized branes in 6D warped 
flux compactiflcations," Phys. Rev. D 76 (2007) 084031 [arXiv:0707.1426 [hep-th]]; 

S. A. Appleby and R. A. Battye, "Regularized braneworlds of arbitrary codimension," Phys. 
Rev. D 76 (2007) 124009 [arXiv:0707.4238 [hep-ph]]; 

C. Bogdanos, A. Kehagias and K. Tamvakis, "Pscudo-3-Branes in a Curved 6D Bulk," Phys. 
Lett. B 656 (2007) 112 [arXiv:0709.0873 [hep-th]]; 

O. Corradini, K. Koyama and G. Tasinato, "Induced gravity on intersecting brane-worlds Part 
I: Maximally symmetric solutions," Phys. Rev. D 77 (2008) 084006 [arXiv:0712.0385 [hep-th]]; 

F. Arroja, T. Kobayashi, K. Koyama and T. Shiromizu, "Low energy effective theory on a 
regularized brane in 6D gauged chiral supergravity," JCAP 0712 (2007) 006 [arXiv:0710.2539 

[hep-th]]; 

O. Corradini, K. Koyama and G. Tasinato, "Induced gravity on intersecting brane-worlds Part 
II: Cosmology," Phys. Rev. D 78 (2008) 124002 [arXiv:0803.1850 [hep-th]]; 

V. Dzhunushaliev, V. Folomeev and M. Minamitsuji, "Thick brane solutions," arXiv:0904.1775 
[gr-qc]. 

[16] K. Lanczos, Phys. Z. 23 (1922) 239-543; Ann. Phys. 74 (1924) 518-540; 
C.W. Misner and D.H. Sharp, Phys. Rev. 136 (1964) 571-576; 
W. Israel, Nuov. Cim. 44B (1966) 1-14; errata Nuov. Cim. 48B 463. 

[17] T. Damour and A. M. Polyakov, "The String Dilaton And A Least Couphng Principle," Nucl. 
Phys. B 423 (1994) 532 [arXiv:hep-th/9401069]; 

P. Brax, C. van de Bruck, A. C. Davis and D. Shaw, "The Dilaton and Modified Gravity," 
Phys. Rev. D 82 (2010) 063519 [arXiv:1005.3735 [astro-ph.CO]]. 

[18] E. Dudas, C. Papineau and V. A. Rubakov, "Flowing to four dimensions," JHEP 0603 (2006) 
085 [arXiv:hep-th/0512276]; 

E. Dudas and C. Papineau, "Dual realizations of dynamical symmetry breaking," JHEP 0611 

(2006) 010 [arXiv:hcp-th/0608054]: 

C. P. Burgess, C. dc Rham, L. van Nicrop, "The Hierarchy Problem and the Self-Localized 
Higgs," .JHEP 0808 (2008) 061. [arXiv:0802.4221 [hep-ph]]. 

[19] S. Weinberg, Gravitation and Cosmology, Wiley 1973. 

[20] C. W. Misner, J. A. Wheeler and K. S. Thorne, Gravitation, W. H. Freeman & Company 1973. 

[21] A. Vilcnkin, "Gravitational field of vacuum domain walls and strings," Phys. Rev. D 23, 

(1981) 852-857 

W. A. Hiscock, "Exact gravitational field of a string," Phys. Rev. D 31, (1985) 3288-3290 



-34- 



[22] A. J. ToUey, C. P. Burgess, C. de Rham and D. Hoover, "Scaling solutions to 6D gauged chiral 
supergravity," New J. Phys 8 (2006) 324 [arXiv:0608.083 [hep-th]]; 

A. J. Tolley, C. P. Burgess, C. do Rham and D. Hoover, "Exact Wave Solutions to 6D Gauged 
Chiral Supergravity," JHEP 0807 (2008) 075 [arXiv:0710.3769 [hep-th]]; 

M. Minamitsuji, "Instability of brane cosmological solutions with flux compactifications," 
Class. Quant. Grav. 25 (2008) 075019 [arXiv:0801.3080 [hep-th]]. 

H. M. Lee and A. Papazoglou, "Codimension-2 brane inflation," Phys. Rev. D 80 (2009) 
043506 [arXiv:0901.4962 [hep-th]]. 

[23] H. M. Lcc and A. Papazoglou, "Scalar mode analysis of the warped Salam-Sezgin model," 

[hep-th/0602208]; 

C. P. Burgess, C. de Rham, D. Hoover, D. Mason and A. J. Tolley, "Kicking the rugby ball: 
Perturbations of 6D gauged chiral supergravity," JCAP 0702 (2007) 009 
[arXiv:hep-th/0610078]; 

S. L. Paramcswaran, S. Randjbar-Dacmi and A. Salvio, "Gauge fields, fermions and mass gaps 
in 6D brane worlds," Nucl. Phys. B 767 (2007) 54 [arXiv:hep-th/0608074]; 

"Stability and Negative Tensions in 6D Brane Worlds," JHEP 0801 (2008) 051 
[arXiv:0706.1893 [hep-th]]; 

"General Perturbations for Braneworld Compactifications and the Six Dimensional Case," 
JHEP 0903 (2009) 136 [arXiv:0902.0375 [hep-th]]. 

[24] W. D. Goldbcrgcr, M. B. Wise, "Renormalization group flows for brane couplings," Phys. Rev. 
D65 (2002) 025011 [arXiv:hep-th/0104170]; 

E. Dudas, C. Papineau, V.A. Rubakov, "Flowing to four dimensions" JHEP 0603 (2006) 085 
[arXiv:hep-th/0512276vl]; 

C. dc Rham, "The Effective Field Theory of Codimension-two Branes," JHEP 0801 (2008) 

060 [arXiv:0707.0884 [hep-th]]. 

F. del Aguila, M. Perez- Victoria, J. Santiago, " Effective description of brane terms in extra 
dimensions," JHEP 0610 (2006) 056 [arXiv:hep-ph/0601222] 

K. A. Milton, S. D. Odintsov, S. Zerbini "Bulk versus brane running couplings" Phys. Rev. 
D65 (2002)065012 [arXiv:hep-th/0110051] 

[25] C. P. Burgess, "Quantum gravity in everyday life: General relativity as an effective fleld 
theory," Living Rev. Rel. 7 (2004) 5 [arXiv:gr-qc/0311082]. 

J. F. Donoghue, "Introduction to the Effective Field Theory Description of Gravity", 
[arXiv:gr-qc/9512024] 

[26] C. P. Burgess, P. Grenier and D. Hoover, "Quintessentially flat scalar potentials," JCAP 0403 

(2004) 008 [arXiv:hep-ph/0308252]. 

[27] Katherine Freese, Joshua A. Frieman and Angela V. Olinto, "Natural inflation with pseudo 
Nambu-Goldstone bosons" Phys. Rev. Lett.65 (1990) 3233—3236 . 

[28] C. M. Will, "The confrontation between general relativity and experiment," Living Rev. Rel. 9 

(2005) 3 [arXiv:gr-qc/0510072]; 

E. G. Adelberger, B. R. Heckel and A. E. Nelson, "Tests of the gravitational inverse-square 
law," Ann. Rev. Nucl. Part. Sci. 53 (2003) 77 [arXiv:hep-ph/0307284]. 



- 35 - 



